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Message  from  the  Chair 


The  record  heat  wave  in  Japan  did  not  dim  the  spirits  or  the  technical  quality  of  the 
ISS AC’90  conference  which  was  held  in  Nihon  University,  Tokyo  (August  20  -  24).  The 
188  participants  come  from  many  different  countries  including  a  good  number  from  the 
USSR.  The  conference  also  enjoyed  strong  local  support  in  terms  of  attendance  and 
industrial  contribution.  In  parallel  with  the  technical  presentations  were  10  demos  on 
various  workstations  through  the  entire  meeting.  It  is  worth  noting  that  ISSAC’90  marks 
several  FIRST’S:  first  time  the  Symposium  ever  comes  to  the  Far  East,  first  use  of  the 
ISS  AC  logo,  and  the  first  issue  of  ISS  AC  Proceedings  to  be  distributed  by  Addison-Wes- 
ley. 


Furthermore,  a  limited  number  of  copies  are  available  at  ACM  member  discount: 

$18.00  for  members  (ACM  order  number  505900,  ACM  ISBN  0-89  791-401-5).  Send 
order  with  payment  to  ACM  Order  Dept.,  P.O.  Box  64145  Baltimore,  MD  21264;  credit 
card  orders  (1-800-342-6626).  Non  members,  please  order  from  Addison-Wesley 
Publishing  Company,  Order  Dept,  Jacob  Way,  Reading,  MA,  01867;  1-800-447-2226 
(A-W  ISBN  0-201-54892-5). 

Now  that  the  ISSAC’90  is  behind  us  we  can  look  forward  to  the  next  meetings. 

ISSAC’91  (Bonn,  Germany!)  has  announced  a  paper  of  deadline  of  Dec.  31,  1990.  Let 
me  encourage  everyone  to  prepare  and  submit  your  contributions  early.  The  proceedings 
will  again  be  distributed  by  Addison-Wesley. 

The  decision  was  made  in  Tokyo  that  University  of  California  at  Berkeley  will  host 
ISSAC’92.  Richard  Fateman  volunteered  as  Local  Arrangements  Chair.  The  organiza¬ 
tion  of  the  conference  and  program  committees  is  underway.  A  primary  focus  of 
ISSAC’92  will  be  the  interface/integration  of  symbolic  and  numeric  computation.  There 
really  should  not  be  artificial  boundaries  between  symbolic  and  numeric  computations. 

Instead  we  should  think  of  mathematical  computation  for  science  and  engineering  as  a 
whole.  This  view  is  gaining  support  among  the  SIGSAM  and  SIGNUM  membership. 

In  fact,  the  two  SIG’s  plan  to  join  forces  at  ISSAC’92.  If  you  are  interested  in  helping 

with  this  conference,  please  feel  free  to  contact  me  or  others  responsible  for  ISSAC’92.  - - - 

an  For 

Preparations  for  the  election  of  new  SIGSAM  officers  are  in  high  gear.  The  *4*1 
nominating  committee  consists  of  R.J.  Fateman  (chair),  A.C.  Hearn,  R.  Jenks,  and  myself.  ^ 

We  are  making  good  progress  in  contacting  people  and  getting  commitments  to  run  for  lced 
office.  Hopefully,  statements  from  the  candidates  can  also  appear  in  the  January  issue  11  orL 
of  the  Bulletin. 
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Table  of  Special  Function  Integrals 


Jean  C.  Piquette 


Naval  Research  Laboratory,  Underwater  Sound  Reference  Detachment, 
P.0.  Box  568337,  Orlando,  Florida  32856-8337 


This  paper  is  the  result  of  a  suggestion  by  B.  F.  Gaviness  [1]  to  publish 
a  table  of  integrals  obtained  using  the  special  function  integration  technique 
that  has  been  described  in  detail  elsewhere  [2] -[8],  since  only  a  much 
abbreviated  version  of  the  table  could  be  accommodated  in  [7] .  This 
integration  technique  permits  symbolic  evaluation  of  integrals  of  the  form 


a) 


f  00  n  «  Ci) 

i=l  'l 


(x)  dx  , 


where  f (x)  is  an  essentially  arbitrary  function  and  IIR  is  a  product  of  any  of 
a  very  broad  class  of  special  functions  including,  but  not  limited  to,  Bessel 
functions,  Legendre  functions,  Hermite  functions,  Laguerre  functions, 

Chebychev  polynomials,  Jacobi  polynomials,  and  Gegenbauer  polynomials.  In  [7] 
is  described  a  method  which  can  be  used  to  enhance  the  on-board  integrator  of 
a  symbolic  software  system  to  include  analytical  symbolic  evaluation  of 
integrals  of  this  type.  The  method  is  currently  being  evaluated  for  possible 

inclusion  in  the  program  Matheroatica  [9] . 

The  table  is  devided  into  two  sections.  In  Sec.  I  are  presented  results 
derivable  directly  from  the  technique  described  in  [7] .  In  Sec.  II  are 
presented  results  that  were  derived  using  the  methods  of  [2] -[8],  but  these 
examples  required  human  intervention  for  their  evaluation.  Current  research 
in  the  area  of  this  integration  technique  falls  into  two  broad  categories  that 
roughly  correspond  to  the  subdivision  of  the  tables.  The  first  area 
corresponds  to  looking  for  heuristics  beyond  those  presented  in  [7]  for 
obtaining  particular  solutions  to  the  ordinary,  linear,  inhomogeneous 
differential  equation  produced  by  the  method.  (Any  particular  solution  of 
this  differential  equation  is  sufficient  to  produce  an  analytical  expression 
for  the  associated  special  function  integral.)  The  second  area  is  associated 
with  looking  for  algorithms  that  can  eliminate  the  need  for  human  intervention 
to  produce  results  of  the  type  presented  in  Sec.  II  of  the  tables.  Current 
efforts  center  around  attempts  to  automatically  incorporate  the  use  of 


generalized  hypergeometric  functions,  F  , 
particular  solution  search  for  these  cases. 


or  Meijer’s  G-function, 


into  the 


In  the  tables,  Z  and  2  denote  Bessel  functions  J  or  Y,  P  denotes  Legendre 
functions  of  either  the  first  or  second  kind,  H  denotes  Hermite  function,  and 
L  denotes  Laguerre  function.  The  parameters  fi  and  v  are  arbitrary  numbers. 
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I.  Results  obtainable  by  fully  automatic  methods. 
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+  ~4  x  Pl/2^x)  tP3/2^x^ 

* •(-  is  - !  x2} 

*  dgfs  *  3x3)  [P^Wl^W 
"  S  fP3/2  M  J 
2  -x2 

"x  H  (x)  fi^xjdx  =  5^  [-^(x)  *  E^+1  (x)  B„(*)] 


f  __•>  e"x  f  Vx)  Vx> 

J  xe  H^(x)  ®j/(x)  ^  ~  2  (l-^+l/)  (l+/i-i/) 

+  V1(X)  Sj/^  d+;-)  +  H ^  ®*/+l^  (l-^) 


Wx>  Wx> 

(1-/B+J/)  (1+^-1/) 


2  2 
}  xVx  H  (x)  B^x)  dx  =  e"x 


2  [-H.W  ^(x)  X (/*+!/) 


(2-fi+v)  (2+p-v) 


+  Vi(x)  Bi/(x) 
-  y*>  *Wx) 


f2^+3i/+2tfx2-2i/x2-tf2x2-i/2x2+2jitoc2) 
2(p-v)  (2-fi+v)  (2+fi-u) 

f2+3tf+l/-2jix2+2i/x2-tf2x2-t/2x2+2tfioc2^ 
2  (/*-*✓)  (2-p+i/)  (2+fi-v) 


B/i+l^  ®*/+l^  (2-/i+l/)  (2+/*-j/) 
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/  e"3x2  x2  H3(x)  dx  =  e"3x2  ^  x  (5+6x2)  H3(x) 


+  l  (l+6x2)  H2/,(x)H(;/q(x) 


2/3 K  J  5/3’ 


~  8  x  H2/3^H5/3^ 


16  H5/3^x^  ) 


(28) 


I  e  X  L^(x)  E^(x)  dx  =  e-x  [l^(x)  E^x)  +  L^x)  E|/+1(x) 

"  Wx>  Vx)] 


(29) 


2  .2 


/  xe'X  V(x)  E,W  dx  =  e-  [-  ^  ^ 

-  (i^(^sr)  V*>  Wx> 

-  Vi(x)  Vx> 

-  Vi(x)  Wx>]  (3°) 


/  e_3X  x  L2/3(x)  d*  =  e3X  {  L2/3(x) 


853x 


16 

*  L5/3(*> 


675x3  +  225x4 


2/3' 
5 


125 

24 


625x 

24 


16 


125 


16 

125x 


27x 

16 


125x 


+  3 

+  3 


125 

24 


24  12 


125x  275x 

- 7T-  + 


16 


16 


75x 

16 


3  1 


125  125x 

24  6 


515x2  +  135x3 


L2/3(x)L5/3(x) 

4 


16 


8 


45x 

16 


L2/3^r^L5/3^x^  }  (31) 
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II.  Results  requiring  human  intervention.  Here,  £.  denotes  integer. 
/  xP  Z  (x)  Z^(x)  dx  =  AqqCx)  Z  (x)  Z„(x)  +  AQ1(x)  Z  (x)  ^^(x) 


where 


*  Aiof*)  VlW  M*>  *  All«  V»1(X)  Wx)  -  (32) 


a  = _ * _ j)3»  +  l?.+ii+yX  p2^ 

A00  2(/t+i/)  u  A11  2(ji+i/)  U  A11 

(-7-3u-3v-2ui/+U2+i/2-4x2)  n. 

2x(ji+v)  "All 

,  (-2-u-v)  (-4-2»v+h2+v2-2x2)  .  x^; 

2xV*)  11  ^  ' 


i  —  ~x  n^A  +  —  ^  —  r\2 » 

01  '  20.2-2)  11  2(,2-,2)  11 

_  (7-3»2-l/2.4x2)  „a  (4.^2-3)I2-|.3-k2.2x2)  a  «^2 

2(/-*2)  11  x(,2-,2)  11  0-2-,2)  ’ 

A  = _ *! _  D3a _ 3J_  D2a  .  I7-|.2-3A4x2)  „a 

10  2(,2V)  11  20*2-,2)  11  2(,2-,2)  11 

(4+flVtf2-3i/2-t/3+2x2l  .  x^2 

,  2  2,  11  ”  ,  2  2v  » 

x(/i  - v  )  (ji  -v  ) 


A  .  •.  £.+3  T  ,  2n 

A11(x)  =  x  L  dQx  , 

n=0 


where 


d.= _ 21^ _ 

4)  (£.+3)  4-8  (£.+3)  3+2  (12-/*2-*/2)  (£.+3)2-8(£.+3)  (4-/42-i/2)  +  [  (2-/<)  2-^2]  [  (2+/*)  2-i/2] 
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{  (3+2n+£.)  4-8  (3+2n+2.)  3+2  (12 -p2-v2)  (3+2n+£)  2 

-8  (3+2n+£.)  (4 -fi2-i/2)  +  [  (2-/0  2-v 2]  [  (2+/i)  2-i/2] }  d 


=  -4(l+2n+£,)  (2n+2.)dQ_1 


n  n  >  0  , 


and  dQ  =  0  if  n  >  n  .  The  parameter  n  is  given  by 

’  0  if  t  =  -1 

-TS--  -  1  if  i.  <  -1  and  even, 

i  ^ 

n  =  ' 

-■"2 ~  if  £.  <  -1  and  odd, 

»  if  a  2  0 

If  IL^O,  dQ  in  the  sum  for  can  be  written  in  the  form 

2(-4)°(fc+l)2,,,1 

an  n  * 

n  g(k, 
k=0 

where  (£.-*-1) 2n-i-l  *s  a  shifted  factorial,  {£-+1) 2n+ ^  =  (£-+1)  (£+2)  . . .  (2n+£.+l)  ,  and 
g  is  defined  by 

g(k, £,/*,!/)  =  (2k+£+3) 4-8 (2k+€.+3) 3+2 (12-/*2-!/2)  (2k+£+3) 2 
-  8  (2k+£+3)  (4-/*2-i/2)  +  [  (2 -fl)  2-u 2]  [  (2+/0  2-i/2]  . 

I  X1-  Z„2(x)  dx  -  *00 (x)  Z„2(x)  <■  2A01(x)  Z^CxjZ^jCx)  *  An(x)  Z2+1(x),(33) 
where 


W*>  -  \  »2*uw  -  DAnW 


*  An(x) 


*01  (»)  =  5  0An(x)  -  An(x)  , 


Alx(x)  =  x  y(x) 
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=  I 


°p_1 

=  1/22. 


b  = 
n 


2(n.l) [x2-(n.l)2]bntl 


Y  —  < 

if  0  £  n  £  — n~  and  2.  ^  3  ,  but  positive  odd  integer  only. 


Also, 


y  =  J:  F  rki 

y  2 2.  4*l[  2  »  2 


-2.-29+1  -L+2V+1  .  -2.+3 

2  »  2  »  1;  2 


x  J 


if  ^negative  odd  integer,  ^integer  but  /  0, 


_  2x*'Z  P  f gi2  , .  jL+3 

y  -  9  9  9*?  19  9  > 

(i+l)[(2.+l)  -4/]  1  ^  J 


H-2i/+3  H+2u+3 


21 

;  -x  J» 


if  2^0  or  even  integer  (positive  or  negative)  but  2.*2l*/negative  odd  integer. 


/  x  P^CxJP^Cx)  dx  =  AqqCx)  P/<(x)Pi/(x) 


*  Aoi  y*)pi/*i(x)  *  aio  Vi(x)pi/(x) 

*  An  ViWpWIw  . 


where 


Aoo<x>  -  lijb  £W  ‘  (M-uffujtTFj  M« 

+  8000  W  Aii<x)  +  «ooi (x)  DAn«  +  8002  W  d2ah(x) 
+  «003(x)  d3ah(x)  +  8qq4 (x)  D4Au(x) 
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*011 M  ” 


-  [2-3/l-l'+3/ix^+Z'x^+3/i^x^ 

2  2  „  2  2  „  2, 

+  1/  x  -3u  -1/  -6x  1 


2(l+/i)  (ji-v)  (1+p+v) 


,  % _ 3x(-l+x  j 

g012'‘X'  ~  (l+/i)  (yi— i/)  (1+/1+I/) 


*013 w  “  2(1+/i)(/i-i/)  (l+^+i/) 


r  %  _  VX[l-ll+V)  I’Z  +  H  +  U) 
g100w  “  (l+i/)  (/*-*/)  (1+/1+1/) 


8imW  = 


[2-/l-3l/+/ix2+3i/x2+/i2x2 

,  2  2  2  „  2  „  2-, 

+  3i/  x  -ji  -3i/  -6x  1 


2(1+1/)  (/i-l/)  (1+/1+I/) 


r  > _ -3x(-l+x  j 

g102^x'  “  (l+i/)  (ji-i/)  (1+/1+I/) 


,  %  _  -f-l+x^r 

g103w  "  2(1+1/)  0*-i/)  (1+/1+J/) 


An(x)  =  E-  b  x* 
P  K 


The  quantities  b  in  the  sum  for  are  constants  given  by 


b  = _ (2g.)  (1+fQ  (l+i/) 

£'"1  (/!+!/+£.+ 1)  Ol+I/-«.+  l)  [  0 1-1/)  2-^2] 


bp  =  -(p+l)(p+2){(p+3)(p+4)bp+4(l-<5p^_3) 
+  2b  0+i/+/x2+i/2-  (p+2)  2]  > 


Q1+1/+ p+2)  (/i+l/-p)  [(/*-l/)2-  (p+1)2] 


(0  <  p  <  £-1) 
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The  prime  on  the  summation  for  signifies  that  p=0,  2,  . ..,  2.-1  when  2.  is 
odd,  and  p=l,  3,  2.-1  when  2.  is  even. 

/  xV^x)]2  d*  =  ^(x)  pyx)]2  *  2A01(x)PI,(x)P1/+1(x) 

*  Au(x)  [Pwl(x)]2  ,  (35) 

where 


acxj(x)  =  H-~¥  °2ahCx)  +  ("1+x2)  {22U)  *  DAhCx) 
00  2(1+J/;  11  (1 +V)Z  11 

+  Lv+jlI  k  M 

(1+i/)  *11 W 


^1^  "  2(1+1/)^  DAllW  "  ^iW 


2.-1 

2 


An(x)  =  l  agn  x-  (2.  =  positive  odd  integer)  , 


2n 


n=0 


and 


a  = _ 2(^l)2 _ 

£_1  (2.-1)  (2.-2)  (2.+3)-2(2L-l)  (2i/+2l/2-3)-4i/(l/+l) 

(2n+2)  {2  [  (2n+l)  (2n+3)  +  (1-2*/-2j/2)  ]  a2n+2+  (2n+4)  (2n+3)  (<52n  £_g-l)  a2n+4> 

^  —  I  »■  ■  ■  I  ■  .1  I  1  I 

2n  2n(2n-l)  (2n+4)-4n(2*/+2*/2-3)-4t/(i/+l) 

2.-1 

0  in  <  ~2~  ,  2.  =  positive  odd  integer 

and  a2n=0  *or  n  ^  (J2.— 1)  /2 .  The  index  n  may  assume  the  value  zero  only  if 

2.  ^  1.  If  2^1,  the  series  for  reduces  to  a  single  term  that  is  given  by 
the  expression  for  a^^  above. 

If,  fb=0,  u= one  half  odd  integer,  but  /  “1/2,  then 


*n  M  -  \  V2°_1 
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where  n  =  *(2z/+l)/2,  and 


bn+2(2n+l)  (2n+2)  (2n+3)-2(2n+l)  [(l+2n)  -2l/(^+l)]bn+1 

-2n[(2l/+l)2-4n2]b  =2(l/+l )25  n  ;  n=0,  1, 

n  n  j  u  / 

*(2i/+n 


and  b  =0  for  n  >  0 

n  Z 


The  preceding  expression  represents  — simultaneous  algebraic  equations 
in  the  unkn0wns  bp  bg,  b*(2l/+l)/2 
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